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Abstract 

In the framework of the Connes-Lott model based on noncommutative geometry, 
the basic features of a gauge theory in the presence of gravity are reviewed, in order 
to show the possible physical relevance of this scheme for inflationary cosmology. 
These models naturally contain at least two scalar fields, interacting with each 
other whenever more than one fermion generation is assumed. In this paper we 
propose to investigate the behaviour of these two fields (one of which represents 
the distance between the copies of a two-sheeted space-time) in the early stages of 
the universe evolution. In particular the simplest abelian model, which preserves 
the main characteristics of more complicate gauge theories, is considered and the 
corresponding inflationary dynamics is studied. We find that a chaotic inflation is 
naturally favoured, leading to a field configuration in which no symmetry breaking 
occurs and the final distance between the two sheets of space-time is smaller the 
greater the number of e-fold in each sheet. 
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1 Introduction 



The early universe is certainly the best laboratory to probe the structure of space and time 
at very short distances, and during this period of the universe, inflation plays a central 
role in solving important problems occurring in standard cosmology. It is usually assumed 
that the dynamics of scalar fields is responsible for inflation, by producing an effective 
cosmological constant term in the Einstein equations. This leads, as well-known, to a 
quasi de Sitter epoch during which the universe scale factor grows almost exponentially. 
The knowledge of quantum field effects, in this era, is particularly relevant because they 
are the origin of the subsequent structure of the present universe. Thus, all models in 
which some scalar fields are present can be used in principle to describe inflation, provided 
that the corresponding dynamics is appropriate. 

In this respect, the recently proposed gauge models inspired by noncommutative geom- 
etry are a relevant example of theories in which scalar fields and corresponding interaction 
potentials emerge naturally from geometry. In this paper we study, for a wide class of 
these models, the predictions on the resulting inflationary scenarios. 

The programme of noncommutative geometry [jl| shifts the stress from space and 
time, traditionally seen as a geometrical ensemble of points, to the (commutative) algebra 
of continuous complex valued functions defined on it, with the perspective of possible 
noncommutative generalizations. One of the most interesting and promising products 
of this programme is the formulation of the electroweak standard model by Connes and 
Lott PI, §], generalized to include Grand Unified Theories in [Q. Loosely speaking, this 
formulation is based on a doubling of space-time, considered as a two sheeted manifold. 
The algebra of continuous functions and the Hilbert space of fields is thus doubled as 
well. The nontrivial feature of the theory is that the Dirac operator, as well as the gauge 
potential (connection), have some nondiagonal elements, which couple the two sheets of 
space-time. These are classical scalar fields: one is related to the component of the metric 
in the discrete direction, and thus to the distance between the two sheets of space-time, 
and the others are Higgs fields, responsible for the breaking of the symmetry. We will 
discuss in more detail the field content of the theory in the first sections. 
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The action is as usual the square of the curvature of the connection. As the Dirac 
operator contains the metric structure of a manifold, the presence of this coupling gives 
a theory of gravity as well ^. By simply postulating that at the level of two-sheeted 
algebra the action is the square of the curvature one obtains the Einstein-Hilbert Yang- 
Mills action containing naturally an Higgs potential. 

Interestingly, at the classical level this scheme generally predicts relations between the 
physical masses of particles, such as the ratio of the top quark mass over the Higgs mass, 
and fixes in some ranges the values of relevant constants appearing in the electroweak 
standard model 0. 

The aims of such a programme are, in the long run, the prediction of physical masses or 
other experimentally testable quantities, and even more the creation of a unique frame- 
work to include gravity. In this initial stage, however, the stress is more on the un- 
derstanding of how the noncommutative structure of space-time might account for the 
already known features of this (exceptionally well working) electroweak model how we 
understand it presently. This understanding, combined with the tools of noncommutative 
geometry, should then provide the framework for the unveiling of deeper symmetries of 
the model, and the possibilities to go further towards the unification with gravity. 

At the moment all studies have been performed at the classical level, and quantiza- 
tion can only be achieved with the usual field theory tools. This quantization, however, 
although it does not respect the relations between the coupling constants 0, under the 
renormalization flow does not change considerably these relations over various orders of 
magnitude of the scale. A fully satisfactory achievement of the programme would be 
a quantization made on the level of the algebra, and therefore intrinsically related to 
noncommutative geometry. This remains one of the main challenges of the programme. 

On more phenomenological grounds, there is also a problem of scales which requires 
further studies. The scalar field which is part of the Riemann connection and is related 
to the distance between the two sheets of space-time, pertaining to the intrinsic structure 
of the latter, has as a natural scale, the Planck length. On the other side, according to 
the Connes and Lott model 0, at low energy this distance must be the inverse of the top 
quark mass, which is several order of magnitude larger than the Planck length, to correctly 
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reproduce the electrowcak standard model. We would like to stress that an important 
feature of our work is to look at the distance as undergoing a dynamical evolution. It 
remains a challenging programme to study the dynamics responsible for the subsequent 
evolution of the distance down to the electroweak scale. 

In this paper we study the scalar dynamics at the Planck energy since we are interested 
in the cosmological implications of the noncommutative scenarios. As we already noticed 
the natural epoch in which the scalar fields, related to the connections in the discrete 
direction, may play a fundamental role is the inflationary one. Therefore, we have studied, 
for the simplest abelian theory, all the relevant characteristic which a good inflationary 
model should satisfy: a quite large e-fold number N, large enough {N ~ 10^) to solve the 
flatness problem of standard cosmology; an efficient reheating stage, which consists in the 
transfer of the enormous energy density stored in the inflaton flelds into entropy density 
for the relativistic degrees of freedom. 

It is worth-while stressing that, although our analysis has been performed for the 
simplest abelian model, the results are general enough to be looked upon as indicating a 
typical behaviour of a noncommutative inspired gauge model. This is because it already 
contains all the key features of more realistic dynamics based on some non-abelian gauge 
symmetry group. 

We cannot say in fact what the full gauge theory of flelds and particles was at epochs 
shortly after Planck time, and therefore we cannot decide the speciflc model of non- 
commutative geometry to be used, namely which Hilbert space and Dirac operator we 
must use. Nevertheless all considered models present at least the flelds mentioned above, 
namely two scalar flelds which are the discrete components of the gauge and gravitational 
connections. 

This paper is organized as follows: in section 2 we review the model of Yang-Mills 
theory plus gravity in noncommutative geometry, and obtain for the abelian case the 
expression for the bosonic lagrangian. Section 3 is devoted to the study of the scalar 
potential and to the related inflationary dynamics. The one-loop effective potential is 
computed in the adiabatic approximation and the results obtained solving the semiclassi- 
cal scalar fleld equations of motion are discussed. Finally we analyze, at one-loop level, the 
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reheating mechanism by introducing the couphng of the scalars to the fermionic degrees 
of freedom. It is shown that this mechanism is efficient enough to make the transition 
from the vacuum energy density to the radiation dominated universe. In section 4 we give 
our conclusions and outlooks. 

2 General relativity and gauge theories in noncom- 
mutative geometry 

In this section the construction of a Einstein- Yang-Mills theory in noncommutative ge- 
ometry is briefly reviewed. The results, are then applied to the abelian case. The reason 
of this choice is that, as far as the scalar fields are concerned, the nonabelian aspect is 
practically irrelevant (apart of the different irreducible representations to which the scalar 
fields belong to) on the resulting inflationary dynamics. Thus, in this scenario, the abelian 
theory not only represents a toy model for which it is interesting to study the cosmological 
implications, but it also contains the essential ingredients of any more complicate, but 
physically more relevant, gauge theory inspired by noncommutative theory. 

2.1 Gravity 

In order to obtain the Einstein-Hilbert action ||^ we flrst introduce the notion of a i^-cycle. 
A K-cycle |], ^, is a triple {A, H, D) where: ^ is a *-algebra unitarily represented on 
the Hilbert space Ti in terms of bounded operators, and D is an unbounded selfadjoint 
operator with compact resolvent such that the commutator with any element of the algebra 
is bounded. This operator is the generalized Dirac operator. The usual concepts of 
Riemannian geometry are now derived by the i^-cycle. The general idea is the following: 
the Riemannian metric is deflned in terms of the i^-cycle as a hermitian structure < . , . > 
on Vl\,{A), which is the space of 1-forms on ^ [|l], Q and thus is the noncommutative analog 
of the cotangent bundle. As far as the notion of connection is concerned, it is useful to 
think of a connection V on a flber bundle as a linear map 



satisfying the Leibniz rule, where E is the space of sections of the flber bundle and A^(A^) 
is the space of 1-forms on the base manifold M.. In the case of a linear connection on 



V : E 



k\M)®E 
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Ai, observing that we are interested in forms rather than vector fields, the fiber bundle is 
taken to be the cotangent bundle. Furthermore, if on Ai there is a metric (the index E 
hereafter denotes the Euclidean version), a Levi-Civita connection is a linear connection 
satisfying the requirements 

V(7^ = , (2.2) 
T{V) = d-PoV = , (2.3) 

where T(V) is the torsion corresponding to the connection V, and P maps the tensor 
product of two 1-forms into their exterior product. The noncommutative notion of a 
Levi-Civita connection will then be given by (|2l|) with E and both replaced by n]j{A) 

V : n\j{A) ^ n\,{A) ® n\,{A) . (2.4) 

Of the two conditions ( p.2|) and (|2.3| ), the former is naturally implemented by requiring 

d < ^,r] >=<V^,r] > - < ^,Vri > , (2.5) 

whereas the latter can be straightforwardly used. The action of V can be extended to the 
whole universal exterior algebra allowing the definition of the curvature -R(V) as 

i?(V) = -V^ . (2.6) 

We will consider gravity on = 7W x Z2 where M. is the Euclidean space-time. On 
Ai the algebra A is taken as the product of two copies of C°° functions on Ai 

A = C°°{M)®C^{M) . (2.7) 

A acts diagonally on the Hilbert space Ti, whose Dirac operator is 

Let us introduce in fl\){A) an orthonormal basis such that 



<e^,e*^>=5^^ , (2.9) 

7" 
7" 



where N, M = 1, . . . ,5; ■y"' = 7'^e° and e^ej^ = gj^^, with a, /i = 1, . . . , 4 Q The covariant 
derivative operator V is characterized by the matrix of 1-forms defined by 

Ve^ = -<,®e*^ . (2.12) 

Analogously the curvature is characterized by the 2-forms R^.j given by 

R{V)e^ = R^^e^' = {dijZ + u^u;^,)®e^' . (2.13) 

The forms can be written as 



(2.14) 



V 750(0^/ 7^<M 
By imposing the conditions (|2.2| ) and (|2.3|) one can see that can be written in 
terms of the Levi-Civita connection cj^j, on A^, of the metric field (p, and of the auxiliary 
fields llSH. The Euclidean Einstein-Hilbert action is 



= {32nG)-' f Tr[<e^(e^)*]-2A / y^Z^d'x . (2.15) 

JM JM 

Eliminating the auxiliary fields and continuing to the Lorentzian signature one obtains 
^E-H= / ('-T7Z77^-2A + ^5/.(^<9''a) d'^x , (2.16) 



Im \ IQ-nG 2 

where IZ stands for the scalar curvature of M., and cr is a real scalar field such that 



(j) = ijexp{—ka) (/i is an arbitrary mass and k = y/AnG = ^jAir/Mpi). 
2.2 Gauge theory 

In this subsection we use the method just developed to add to the gravity action the 
Yang-Mills action for a nonabelian gauge group p, For simplicity let us consider the 
case SUl^ti) (g) SUR{m) ® UsiX) with L and R standing for lefthanded and righthanded 
chirality, respectively. In this case the space will have the same structure as in last section, 
J\A. = M. ^ Z2 = M.L U M.R1 with M. being a smooth spin manifold. 



*Reniind that the Euchdean Dirac matrices 7^ and 75 obey to the foUowing algebra: 7^^^ = — 7^^, 
{7a., 7.^} = -25u,., 75 =75- 



We will consider the i^-cycle for which Ti is the Hilbert space made of fermions and 
Tp^ belonging to n-plet and m-plet of the gauge groups Uiin) and UR^m), respectively. 
Hence a vector ip oiTi will have the form 

*=(*|) , (2.17) 

h = l,..,ng standing for the family number, with Ug the number of generations. Thus 
the Hilbert space results to be ?i = L^{Sl) © L^^Sr) with S'i(/j) denoting a bundle on 
■Ml{R) of multiplets spinors under the gauge group UL{n){UR{m)). The *-algebra, A, of 
bounded operators acting on 7i, is taken as 

A = C^{M)®{UL{n)®UR{m)) . (2.18) 

An arbitrary element of A, a, results to be a (n + m) x (ra + m) matrix 

/ \ 

a« ' 



(2.19) 



where and are two unitary n x n and m x m matrices of functions, respectively. 
The Euclidean Dirac operator, D, is given by 

with M an n X m mass matrix, and K a matrix responsible for a family number mixing. 

In this formalism the gauge potentials ^4^*^^^^ defined on the manifold M.l{r), are 
contained in the 1-form 

p^E^AoM^iJl,, . (2.21) 

where in ( |2.21| ) we have used the following definitions 

^L(R) ^ ^M;(^)^^/^^af^)9^6f^) , (2.22) 

j 

$ = (j)J2af{Mbf -bfM) . (2.23) 
j 

By imposing the condition 

^af^)6f^) = l , (2.24) 
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it follows that Aj^^^^ defined in ( |2.22| ) are antihermitian n x n (m x m) matrices, which 
are connected to the generators of Uiin) (f/i?(m)). Thus in order to obtain the gauge 
group SUl^u) ® SUji^m) ®Ub{1) one has to require a graded trace condition on p, namely 
Tr{Aj^) = Tr(A^). We can express Aj^^^^ in terms of the generators"'" of SUL{n) (SUR^m)) 

r/ (T^^) (p = 1, .., ^2 - 1; g = 1, .., m2 - 1) 



7 ? / 777 

Ai.-^^,,r^^W'^^-^^,J-^^KB, . (2.25) 



A'^—gnT^-'Wf-^gs. ," , 1„. B„ , (2.26) 

^ V2 v2 y fnln + m) 

where W^^ (Wli^) are the corresponding hermitian gauge potentials and is the Ub{1) 
gauge field. The Euclidean Yang-Mills action, as well-known, is given by the following 
expression 

ly-M = ^ J d'^^ Trie') , (2.27) 

where the 2-form 6 = dp + p'^, dp = bj], and A/" is a normalization constant 

to be determined. From the above considerations we get the following expressions for the 
matrix elements of 6 



e 



LL 
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RR 



\r''Gl + VlAl-Y.a'^Uh'^ + {^^^-<P'Y.<^^MM^b^^ KK^ , 
l^'Gl + V^AJ - E «f + (^^^ - 0' E afM^Mbf]^ K^K , 



Olr = 75 [-^$ - k^^cr + $4^ - 4^$] K , 

Orl = 75 - A;$^^(T - 4^$"^ + $"^4^1 ;^"f . (2.28) 



Note that, for sake of simplicity, in Eqs. ( p.28| ) we have shifted the field $ with respect 



to the definition ( |2.23| ), and we have denoted with 7^^^ = [7^,7'^]/2. Moreover, given the 
covariant derivative ( p.l2|) , we define the gauge covariant derivative as "^liR) = ^^+^L(i?) 



and the Yang-Mills tensor G^f,^) = \/LiR)ji^LiR) _ yL{R) jij.^(R) ^ convenient to 

introduce the following auxiliary fields 

Xhr, - -V-,(^)AJ(^)+Earn&r ' (2-29) 



tPor the generators we assume the normahzation conditions Tr{TpT,^) = 25ppi and TriT^T^) 



25qqi. 
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(2.30) 
(2.31) 



Thus, substituting in { ^TI^ we get 



(^RR 

Qlr 

(^RL 



2 

75 



^$ - k<^^a + $4^ - 4^$] K , 



75 



(2.32) 



The Euchdean bosonic lagrangian density corresponding to the Yang-Mills contribution 
( |2.27D is then given by 



B 



Af 



2^" 



KK^ (^$$t _ ^^Yl) - Xl\ + -Tr [K^ K ($^<l> - ^^Yr) - Xr 



^/lu^L I /^fj-iy /^R 



R. ^ fiv 



(2.33) 



The continuation to the Lorentzian signature is easily obtained by substituting C 



B 



and xa =^ it. Replacing in ( p.33| ) the fields A^^^^ with the physical gauge potentials WJ^^ 
W!l^ and B^. we get 



C 



B 



-Tr 



-Tr 



K^K (^^^^ - (t)'^YR^ - Xr 



KK^ ($<|.t _ ^Vl) - Xl 
+ Tr (k^K) Tr [(1)^$)"^ (P^<l>) + kd^,{<^^^) d^'a + k'^^^^d^ad^'a 
- \n, {glW^lWf,i + glWi;:^W^^^ + glB,^B^^)^ . (2.34) 

where VF^^*-^-* and B^^ are the Yang-Mills tensors for left (right) and Ub{^) gauge fields, 
respectively, and 



9l 



9 b 



n — m 



V2 



V2- 



nm{n + m) 



(2.35) 



Note that, in order to have in Eq.( 2.34 ) the usual kinetic terms for the gauge degrees 
of freedom it is necessary to put Ql = Qr = Qb- To obtain the final expression for 
the bosonic action we need to eliminate from ( |2.34|) the contribution coming from the 
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auxiliary fields ^L(i?,) and ^^(/j), and this is performed by varying the lagrangian with 
respect to them and substituting the resulting constraints in ( ^^.34 ). It is immediately 
clear from this procedure that the presence, in the final expression of Cbi of potential 
terms for the $ field is strongly depending on the form of the M matrix. To be more 
precise, it is connected to the possibility, for the particular choice made for M, of not 
having independent auxiliary fields. 



2.3 A simple abelian model 

As clear from the previous sections, the scalar fields content of a noncommutative Einstein 
-Yang-Mills theory is essentially independent of the particular gauge group chosen. Any 
models will have at least a real a field coming from gravity and one or more $ Higgs 
fields. This simple consideration suggests to consider the simplest gauge theory, that is 
n = m = 1, in order to study the resulting cosmological implications. The point we wish 
to make is that the corresponding results will reliably sketch the inflationary scenario of 
any more complicate but more realistic GUT model. 

According to Eq. ( p.20|) , the new Dirac operator in the Euclidean space results to be 



D 



(2.36) 



Consequently, following the previous analysis we get the expression for the matrix elements 
of 6, which read 



Oll = --gBr''B^,-XL+[\^^-iJi^e^^{-2ka))KK^ , 

Olr = -i,m + k<^^^) , 
Orl = -75 + A;$*^a) , 



(2.37) 



where $ is an usual complex field. Thus we get the bosonic lagrangian in Lorentzian 
space-time which is 



-9 



1 



Tr{K^Kf - 



{Tr{K^K)f 



■Tin 



1$^ -/i^ exp(-2A;(T) 



(2.38) 
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From ( |2.38| ) it is clear that only for K'^ K ^ (3 1^^ the potential term for the Higgs field $ 
is preserved. With M = nggB^ we obtain the usual lagrangian. Furthermore, redefining 



1/2 



^ = if qb ng/Tr{K'^ K) , and finally fi = m qb ng/Tr{K^ K) 



1/2 



we get 



B 



4 



X 
4! 



1^1 



m exp {—2ka) — Vq 



+ 



(2.39) 



In ( |2.39| ) the kinetic term of the cr-field and the constant cosmological term appearing in 
( p.l6|) , Vq, are also included, and the definition 



\ = 6gl 



n. 



,Tr{K^Kf 



> 



(2.40) 



is used. From Eq. ( |2.39D one easily gets the tree-level potential for the cp and a fields 



V{a,^) = Vo + 



4! 



vr? exp (—2ka) 



(2.41) 



This potential has its minimum for = exp (— 2/c(t), which is a condition which 
does not fix both \(p\ and a uniquely. The situation changes if one assumes a to be still a 
classical field, and consider the one-loop corrections to V{a, (f) coming from the quantum 
fluctuations of 0]. This assumption is consistent with the physical scenario we want 
to describe, since the inflationary epoch is assumed to be described by Einstein gravity 
(classical theory). In our case this means that the classical Einstein theory can be applied 
to the doubled manifold, so a, which is the metric component in the discrete dimension, 
has to be a pure classical field with no quantum fluctuations. On the contrary, ip is an 
usual scalar field and thus for it we will consider quantum fluctuations. 

3 Inflation from scalar dynamics 

To study the inflation induced in the early universe by the dynamics of a and (p, one 
needs to solve the equations of motion for these fields starting from reasonable initial 
conditions. In this section we will discuss the evolution of this coupled system showing 
how this evolution corresponds, for a natural choice of initial conditions, to an inflationary 
phase. 
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3.1 The effective potential 

The Higgs field contrary to a which we will treat as a pure classical field, can be 
expanded around its background configuration as 



(3.42) 



where the background configuration Lpd = (v?) can always be chosen real using the U{1) 
gauge symmetry of the theory. With 5ipi we denote the independent quantum fiuctuations, 
which satisfy {5ipi) = 0. By substituting ( p.42| ) in ( p.39| ), we get the following equations 
of motion for ipc, o", and 6(fi 



2m^ exp {-2k(T) + (3 Sipl + Sifl) 







(3.43) 



1 + (ifl + {5^1 + 8^1))] Da + 2fcVc<9^<^c<9''fT + kd^^cd^(pc + k^^\J(p, 



+— m^exp (~2ka) 



2m^exp {-2ka) + {5ip\ + 5ipl) 



(3.44) 



+ ^ 3(^^ - 2m^ exp {-2k(T) 5ipi + ^ 5(pl + 5ipi5ipl 



X 

+ 24^^ 



3 (5^1 - {5^1)) + - {5^1) 







(3.45) 



n'^V^2 + ^ '^c ~ 2m^ exp {-2ka) 5lp2 + ^ 5v52 + 5v^?5v'2 + 2Lpc5(fi5Lp2 = . (3.46) 



The quantities {Sip"^) in Eqs. ( |3.43|) and (|3.44| ) account for the one-loop corrections to 



the effective potential. Note that in the quantum fluctuation equations of motion ( |3.45| ) 
and (|3.46|) the equation of motion for has been used, and according to the adiabatic 



approximation ||I0|, we have neglected all terms proportional to derivatives of the back- 
ground flelds. 

In order to linearize the last two equations with respect to 6(pi, one can estimate the 
self-interaction terms of these flelds in the mean-field approximation. According to this 
ansatz, the quadratic and cubic terms of 6ipi result to be 6ipf = 3{S(p'f)Sipi, 6ipj = {Sipf) 
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and 5ipi5ip2 = = respectively. In particular the vanishing of the 6ipi6ip2 corre- 

lator can be easily understood from the explicit expression of the interaction term of the 
fluctuations; from (|2.41| ) 

A 



96 



(3.47) 



one get the conclusion that, at all orders in perturbation theory, there are no contributions 
to graphs with an odd number of 6ip2 external lines so, in particular, {6ipi6(p2) = 0. We 
also notice that all expectation values of cubic terms, like (5v?f), do not contribute to 
( p.45|) and ( p.46| ). By taking the expectation value of these equations, in fact, and using 
the condition (Sipi) = 0, it follows that the sum of all these contributions identically 
vanishes. 

The linearized expressions for Eqs. ( p.45| ) and (|3.46|) take then the form 
A 



24 



3(^^ - 2m^exp {-2ka) + (3 5(pi + 5ipi) 5(pi = 



(3.48) 



n\Sv^2 + ^ [pI - exp {-2ka) + {5^\ + 3 5^1)] 5^2 = Q ■ (3.49) 



By following |T0|, at zero-th order in the adiabatic expansion, we get the gap equations 
for {S^f) 



X 



log 



A (3^1 - exp {-2ka) + ^5v\) + {5^1) 



A [3(/?2 _ 2rn2 exp {-2ka) + 3(5^?) + {6^1)] 



96g2 



(3.50) 



{5^1) 



X 



log 



A [^l - 2m2 exp {~2ka) + {6^1) + ?>{5lpI)] 



96g2 



+ 1 



(3.51) 



Note that for simplicity, we have neglected the contribution of the coupling of ip with 
the scalar curvature i?, because it rapidly becomes negligible. Thus the Eqs. (|3.50| ) and 
( p.51|) are just the flat space corrections a la Coleman- Weinberg 0. 

The parameter Q is the ultraviolet cut-off which, representing the energy upper bound 
for the validity of our model, will be assumed in the following of the order of the Planck 
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mass. These two expressions can be solved by iteration. In the following we will only 
consider the zero-th order approximation, which already gives results correct at order A 

1 



X 



log 



A [Sifil - 2m2 exp {-2ka) + QV(27r2)] 



{5^1) 



96Q2 



(3.52) 



log 



A l^l - 2m2 exp {-2ka) + QV(27r2)] 



96g2 



+ 1 



(3.53) 



From Eqs. (|H) and (HD one reads the square of the effective masses for the quantum 
fluctuations 



24 
X_ 
24 



3(/?^ - 2m^ exp (-2A;a) + 



(fl — 2m^ exp {—2ka) + 



01 

27r2 
27r2 



(3.54) 
(3.55) 



As it always occurs in scalar field theories, these masses get contributions from the highest 
scale involved in the model, which in this case is the Planck mass. 

By virtue of Eqs. (|3.52| ) and ( p.53| ) one can determine the quantum fluctuation con- 
tributions to (|]4|) and ( ^^441) 



5ipc 







(3.56) 



□ a + 2k'^Lpcd^Lpcd^a + kd^Lpcd^ipc + kipcd'^c + 



(3.57) 



where the effective potential Veff{a,(pc), in terms of the definitions ( p.54| ) and (|3.55| ), 
results to be 



A r 



Vo + — 
" 96 



647r2 



(fl — 2rn? exp [—2ka) 



2 , AQ2 



+ 



967r2 



(p\ — rn? exp [—2ka) 



log 



4Q2 



log 



4Q2 



1 

+ 2 



(3.58) 



Notice that in the factor in front of Qcr in Eq. ( |3.57| ) we have only considered the zero-th 
order term in the A expansion; this is to be consistent with the fact that, according to 
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the adiabatic expansion, we have neglected in Eqs. ( p.45| ) and ( p.46| ) the terms order A 



which are proportional to the derivative of the a field. 

The quadratic and logarithmic divergences in V^fj when ^ oo are usually re- 
moved by the renormalization procedure, namely by a redefinition of the mass and the 



coupling constant. This is always possible for a potential of the form ( |2.41| ) if the term 
rn? exp(— 2A;cr) is just a constant. In our case, however, the tree level expectation value of 
the field depends on the dynamical field a . If one assumes that the latter evolves quite 
rapidly with respect to <y9, once a has reached a stationary point cr^, the dynamics of the 
(yj becomes the one governed by the usual mexican hat potential 

V{a,ip) = Vo + ^[\ip\^-m^exp{-2ka^)f , (3.59) 

which can be renormalized in order to take into account the effects of quantum fluctu- 
ations. In the general case we are dealing with, where the simultaneous evolution of (p 
and a is considered, our model is not renormalizable, i.e. both divergences cannot simply 
absorbed in the parameters of the lagrangian density. This fact is not surprising, since the 
a field is related to gravity in the discrete direction, and still we are far from a renormal- 
izable theory where gravitational fields are considered as dynamical degrees of freedom. 
In the following we will look at ( |3.58| ) as an effective potential and we will fix the cut-off. 



as already mentioned, at the order of the Planck mass scale, which is the natural one at 
which the field dynamics takes place. We will see, by numerically solving the one-loop 
equations of motions, that in fact the evolution of a is quite fast due to the exponential 
nature of its potential. However, the largest contribution to the e-fold number just comes 
from the time interval during which a reaches its stationary point, so the a dynamics, as 
far as the infiationary dynamics is concerned, is particularly relevant. 



The study of the minima of the effective potential (|3.58|) can be analytically performed 
in a simple way if we consider only terms up to the first order in A: in this case the only 
minimum is located at the point ipc,m = 0, (Xm = (1/2A;) log(87r^m^/(5^). By requiring 
that the effective cosmological constant vanishes at the absolute minimum we therefore get 
Vq = A(5^/(15367r^). We have numerically checked that the shift in the true minimum due 
to the order A^ logarithmic terms is of few percents for A ~ .1, which will be the largest 
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value we will use in the following in order to have a reliable perturbative treatment of the 
scalar dynamics. 

In Figure 1. the effective potential for a choice of parameters, namely Q = 10 Mpi, 
m = 1 Mpi and A = 0.1, is shown. 

3.2 The equations of motion 

Since we are interested in the cosmological implications of the scalar dynamics, and in 
particular in the possibility to get a reasonable inflationary stage, we choose a spatially 
fiat Friedmann-Robertson- Walker (FRW) universe, whose metric has the simple form 



ds'^ = dt^ - a^{t) dx^ 



(3.60) 



and assume homogeneous field configurations, namely a = cr(t), ipc = ^Pc{t)- 

Under this assumption, the equations of motion Eqs. (|3.56| )- (|3.57| ), once put in normal 

form, read 



a + 3Ha + 1 + 



eg' 

47]-2 



2 , ^Veff ^V^ff 



5a 







(3.61) 



OCT Oifc 

-P^^a' + ^ = 



where H is the Hubble parameter, defined as 

~ -2 



H 



3 1^ + ^^^^'=^+^ 



47r2 



1/2 



(3.62) 



(3.63) 



We will fix the initial conditions for a and <^c by requiring that the corresponding initial 
value for the energy density is of the order of Mpi, with vanishing kinetic terms. Moreover, 
we also require that the initial distance between the two sheets, ^5(0), is equal to the 
Planck length, to be consistent with the idea that our model is able to describe the 
evolution of the universe just after the quantum epoch for gravity. Since all dynamics 
only depends on the combination rf^^ = m exp(— A;o"), and m can be reexpressed in terms 
of ^5(0), we get ^5 ^ = d^^{0) exp[— A;(cr — cr(0))]. This allows, by shifting a ^ a + cr(0), 
to take (t(0) =0 as initial condition without loss of generality, since all the other terms 
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in the equation of motion are invariant under this transformation. Consequently, the 
requirement on the initial energy density fixes, up to an inessential sign, (pc{0). 

For the free parameters of Veff, A and Q, we have considered the following ranges: 
A G [0.01,0.1], which ensures the reliability of the perturbative approach, and Q > Mpi. 
The latter choice is, of course, based on the reasonable assumption that in the initial stage 
the only typical mass scale is Mpi. 

The set of coupled equations ( |3.61|) - (|3.63| ) has been numerically solved. We will here 
sketch the main features of the a — (pc motion: a more detailed analysis is reported in the 
next section, where also the reheating epoch is considered. 

The dynamics is characterized by two main behaviours: in a first phase the scalar fields 
evolves towards the minimum a = (Tm,f ~ (see section 3.1). During this rolling down, 
as in the chaotic scenarios, the scale factor rapidly grows, faster the larger is the value of 
the cut-off parameter Q. It is interesting to note that, since the stationary point for ip 
differs from zero for terms of the order A^, the U{1) symmetry is broken quite weakly. In 
particular, if the choice for Q satisfies the condition Q < 2ttMpi it is easy to see that the (pc 
field undergoes an inverse phase transition. At t = in fact, having chosen ^5(0) = Mpi, 
^liVc = 0) < 0, showing that we are in a broken symmetry phase: in these conditions, 
in fact, the minimum at fixed a = a{0) = is located at (pl^ = 2M|,; - Q^/2tt^. When 
d^, due to the a evolution, reaches the critical value d^^ = 2-n/Q > d^^O) the symmetry 
is restored, up to terms order A^ due to the logarithmic terms in Veff, since the true 
minimum becomes ipc = 0. Differently, if Q > 2'kMpi all the evolution take place in 
the unbroken phase, since in this case d^ always remains above the critical value. 

Once the rolling down stops, the value of a remains practically fixed at dm, while 
starts oscillating around the equilibrium configuration. This phase is extremely important 
for reheating since, as well known, during these oscillations the scalar fields may transfer 
their energy density to the incoherent gas of 5(pi quanta and, more generally, to other 
particles via decay processes. 
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3.3 The Reheating 

It is well known that the necessary last stage of the inflationary dynamics is the so-called 
reheating era. During this phase the energy density stored in the scalar sector is released to 
the relativistic degrees of freedom via decay processes. This epoch is at least as important 
as the truly inflationary one, since it guarantees the end of the exponential-like growth 
of the scale factor and the reheating of the universe, which will then continue to expand 
as in a radiation dominated era. It is usually stated that the infiaton field configuration, 
once it has reached its minimum, can be regarded as a condensate of the corresponding 
massive quanta, which would decay in other light particles. Such a coupling is assumed 
to give rise to an effective term of the form r<i>, where F is the decay rate, which is put hy 
hand into the infiaton equation of motion. It has been shown however that this procedure 



is unsatisfactory |12]. A better approximation consists in deducing the friction terms 



from the one-loop effects of quantum fluctuations at first order in the adiabatic expansion. 



as done in |10| and |12| . This procedure would imply the appearance of terms of the form 
TH(^ in the infiaton equation of motion and, moreover, that the mass of the $ quanta 
should be considered as time dependent. We will restrict ourselves in the following to 
this approach; this will give a correct qualitative description for small times, but, as any 
other perturbative computation, it would fail for large times. This is due, in particular, to 
resonance phenomena which enhance the order % term amplitude over the tree amplitude, 
so one should make use of a Hartree or mean field technique similar to the one discussed 



in the previous section |jT2|. This would lead to a set of gap equations to be solved at 
higher order in the adiabatic expansion. However it is outside the aim of this paper a 
detailed analysis of the reheating of the universe: the final stage of the inflationary era, 
of course, strongly depends on the particular model one assumes, as well on the way the 
infiaton field is coupled to the relativistic degrees of freedom. What we would like to show 
is simply that even the simplest abelian model may give rise to an efficient transfer of the 
energy density from the coupled system of scalar field configurations a — (fc radiation. 

To this end we will describe the reheating phase as due to decay processes of the ip 
field into fermionic degrees of freedom (see eq. ( p. 17] )), due to the interaction term 



Tr[^{D + p)^], where D is defined by (|]3|) and p is the abelian version of the gauge 
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potential matrix ( |2.21| ). These processes lead to an energy transfer from the coherent 
classical configuration, whose evolution drives the inflationary expansion of the universe, 
to the incoherent gas of quanta: if the latter starts to dominate the energy density of 
the universe, the inflation stops, since their equation of state relating the energy density 
pR and the pressure pr satisfies the condition pR + 3pR = ApR> 0. 
The interaction term can be written in the usual Yukawa form 

ll^Yuk = -f a{tYip ipLi'B^ + h.c. . (3.64) 



Using the background splitting (|3.42|) we have 



implying that the fermions acquire an effective mass m^, = fipc/V2. For the decay of 
the 6ipi quanta into a fermion-antifermion pair to be possible the scalar particle masses 
should be greater than twice the fermion mass. The corresponding decay probabilities in 
the 6(pi rest frame can be simply evaluated in the Born approximation and result to be 

- f 

Ti = T{6^i^'^^) =nf—nie{n,-2m^) i=l,2 , (3.66) 

167r 

where ~ 10^ is the number of independent decay channels. In the general case of a 
non-abelian gauge theory, one should also consider interaction terms of the scalar field 
$ with the gauge bosons, see ( |2.34| ). This would increase the number of possible decay 
channels for $ and the efficiency of the entropy release to the relativistic degrees of 
freedom. However the analysis for the simple abelian case will be general enough to show 
how the model may also quite naturally account for the end of inflation. 

The (pc and a equations of motion modify due to the presence of the interaction terms 
with the fermionic field: the expectation value of the fluctuation field squared (Sipf) should 
in fact include the interaction with the \1' degrees of freedom. Following ||T0| , to which for 
brevity we refer for the detailed calculation, one finds 



rr 



Ti (3.67) 



2567r 

where with {5ipf)Q we have denoted the corresponding result in absence of the interaction 
term ( p.65|) , which has been quoted in (|3.52| ) and (|3.53|) . Using these expressions we can 
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compute the equations of motion for a and ipc] as before, inserting ( |3.67| ) into ( p. 431 ) and 
( |3.44 ) and passing to normal form 



a 



3H& 



5o 



- k<^c- 



6(pc 
, 



(3.68) 



+ 3ii - fcv^c I 1 + 



k{ipc + kifcdf + 



5V, 



- kLp, 



eff 



+7,^ - kifclv, - k'^'^ccr'^ + 



5V, 



eff 



+ 7^. = 



With and we have denoted the following quantities 



la 



(3ri + r2)i/ + 3[|ri 



61447r 
Xkim? exp(— 2A;cr) 



VL2 



(3.69) 



(3.70) 
(3.71) 



30727r 

which represent the dissipation terms due to fermion-antifermion pairs production. 

The complete system of differential equations describing the evolution in the inflation- 
ary phase we are interested in, includes also the equation ruling the evolution in time of the 
energy density of the relativistic fluid of ^ field quanta, which are produced via the decay 
term. This equation is provided by the covariant conservation of the energy-momentum 
tensor which reads 



•4>c) 1 



where 



0^ 



1 



Pa,^, = TT + k^pc^cO- + - 



P. 



l + k'hl 



0^ 
QY 

47r2 



~ ^eff{cr, <-Pc 



(3.72) 



(3.73) 



(3.74) 



represent the pressure and the energy density of the total system of a — ipc fields, respec- 
tively. 

Equation ( |3.72|) , by using the definitions of p^^ip^ and Pa,ipc and the equations of motion 
( p. 681 ), ( p. 691 ), can be cast in the form 



PR + AHpji = 7^cr + 7<^^0c 



(3.75) 
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Finally, for the Hubble parameter H we get the expression 

H ^ {^^^ = [p^,^^ + pnf' . (3.76) 

The system of equations ( p.68|) , (|3.69| ), ( |3.75|) and ( p.76| ) allows to compute the e-fold 
number of the inflationary stage, defined as = log [a(t/)/a(0)], where tj is the 
value of time at which the inflation ends. This time can be fixed as the one at which 
the radiation energy pnitf) eventually equals the energy density Pa,ip^(tf) stored in the 
background configuration (pc — a. In terms of the Hubble parameter ( p.76|) , the expression 
for takes the simple form 

N= H{a{t),^{t),pn{t)) dt . (3.77) 
Jo 

The reheating temperature, namely the temperature at the beginning of radiation domi- 
nated era, can then be found by expressing the radiation energy density pnitf) in terms of 
the temperature. If one supposes in fact that a complete thermalization already occurred 
at due to fast gauge bosons mediated interactions among the fermions, we have 

PRitf)=nf^T'^H . (3.78) 

In Table 1 the e-fold and the reheating temperature are computed for some values of 
Q and A, with / = 1, while in Figure 2 we report the trajectories in the a — (pc plane for 
the same choices of the free parameters. 

As already mentioned in the previous section the (T—(fc motion is initially characterized 
by a rolling down toward the minimum. During this phase, shown in Figure 2, the 
radiation energy density remains very small while the e-fold A^ rapidly increases, reaching 
a value very close to the final one. From Table 1, one can see that its value depends quite 
strongly on the choice of Q and A: A^ grows if Q increases and/or A is small. In particular, 
with Q = 5^10 Mpi N is of the order 10^, which is the value needed to solve the flatness 
problem of standard cosmology. 

When the system reaches the absolute minimum of K//, (pc starts oscillating around 
fc,m- For all cases reported in Table 1, for which the fermion coupling constant has 
been chosen to be quite large, / = 1, the reheating temperature, T^h, evaluated when 
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Pacific. = PRi is of the order of 10^^ -h 10^^ GeV . Interestingly this value is the one typical 
of GUT theory energy scales. If the value of / is decreased, the reheating, in the one-loop 
approximation we have considered, becomes less efficient. In fact, already for / ~ .1 
the energy injection from the a — (pc system does not compensate the dilution of the 
radiation energy due to the universe expansion. Therefore, the value of the ratio PrI Pa,^^^ 
reaches a maximum smaller than the unity, and then starts to decrease again. In these 
conditions the universe would continue to evolve in a Po-,(/>c matter dominated regime. This 
also happens for one of the cases reported in Table 1, namely for the choice Q = \ and 
A = 0.1. The corresponding value for the e-fold number, which is however unsatisfactorily 
small, should be understood as the one at the beginning of the matter dominated era. 

4 Conclusions and remarks 

This paper represents a first step towards a systematic study of the cosmological implica- 
tions of gauge theory and gravitation models inspired by noncommutative geometry. We 
have investigated this issue following 0, ||, |^, where gravitation and gauge theories are 
combined in the framework of noncommutative geometry. By virtue of the two-sheeted 
nature of the manifold, different scalar fields appear due to gravity and gauge connections. 
The scalar dynamics naturally provides a chaotic inflationary scenario in which, just after 
the Planck epoch, an initial field configuration (cr(0), V5c(0)) originates from quantum fiuc- 
tuations with a corresponding energy density of the order of Mp;. The evolution from the 
above initial condition is then fixed by the theory. A remarkable aspect of this scenario is 
that the link between the scalar fields, the distance along the discrete dimension, and the 
role played by the gravity in the doubled manifold A4, naturally fixes the initial values 
for a and <^c- This contrasts with the arbitrariness of the usual models. 

In this first analysis we have chosen an abelian gauge theory coupled with gravity, 
which represents the simplest model capturing all the essential ingredients of the dynamics. 
As is well-known, in noncommutative geometry to have a nontrivial scalar interaction 
(scalar potential) more than one fermion generation is needed. 

We have considered the one-loop corrections to the classical potential due to the 



5(pi quantum fiuctuations introducing a cut-off regularization; differently than in |T0 
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a renormalization procedure cannot be applied since the Lp^ vacuum expectation value 
explicitly depends on cr, as explained in section 3.1. Since from our numerical simulation 
it comes out that this field reaches its stationary point cr^ quite soon, a similar procedure 
can be applied in fact for the subsequent evolution in which only the ipc field is still varying 
while a remains frozen at a„i. To study the inflationary phase, however, since the main 
contribution to the e-fold number comes from the epoch in which a is rapidly evolving, 
we were forced to use the one-loop potential explicitly depending on the cut-off parameter 
Q, for which we have considered values of the order of the Planck mass. 

The results show that, for reasonable values of the parameters, the value of N is of 
the right order of magnitude (A^ ~ 100) and the distance of the two sheets of the 
manifold M. remains of the order of the Planck length. It is interesting to note that it is 
the presence of the discrete dimension which is responsible for the inflationary expansion; 
this can be understood by realising that in the model it is cig ^ which play the role of the 
effective cosmological constant. It is possible, with different choices of the parameters, to 
have signiflcantly growing as well, but this can only be done at the price of a much 
smaller e-fold. 

The reheating phase has been studied by considering the contribution at one loop 
to the {5ipj) of the interactions with fermion pairs. This introduces a dissipative term 
into the and a equations of motion which describe, at this level of approximation, 
the energy release of the scalar flelds to the relativistic degrees of freedom. We found 
that, provided the coupling constant is large enough, / ~ 1, a complete and efficient 
reheating is obtained and that the corresponding reheating temperature is of the order 
of 10^^ IQ^^GeV, well below the initial energy scale, which is set by the Planck scale. 
As we already pointed out, the approximation we have used is unsuitable for a detailed 
analysis of the reheating phenomenon. We can nevertheless conclude that there are strong 
indications of the presence of an efficient reheating stage. 

There are at least two subjects which would be interesting to study further. First 
of all a more realistic model, based on a nonabelian gauge symmetry group, should be 
analyzed in more details, to see its implications for the infiationary scenarios. Incidentally 
the presence of couplings of the scalar fields with the gauge bosons should increase, in 



24 



this case, the efficiency of the reheating mechanism. 

It has been shown that in the model we have considered, the distance at the end of 
the inflationary era remains of the order of the Planck length; it remains an open question, 
however, if its subsequent evolution is non-trivial, namely that when the universe continues 
to expand after infiation, its value gradually increases, up to the inverse top quark mass, 
which is the value needed in order to reproduce the electroweak standard model. It seems 
that the idea of looking at as undergoing a dynamical evolution is quite important, 
since this could produce the onset of several phase transitions which took place during 
the expansion of the universe. 
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Table ] 




Q 

{Mpi) 


A 


N 


Trh 
(10^ GeV) 


10 


0.1 


130 


5. 


5 


0.01 


110 


.6 


5 


0.1 


40 


2. 


1 


0.01 


20 


.1 


1 


0.1 


10 


no reheat. 
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Figure Captions 

Figure 1. The effective potential Veff{a, (pc) is plotted for A = .1 and Q — 10 Mpi and m — 
1 Mpi. 

Figure 2. The trajectories in the a — (pc plane corresponding to the values of the parameters 
reported in Tabic 1, are shown. The labels refer to the order in which they are 
considered in the Table, from top to bottom. 



Table Caption 

Table 1. The e-fold number, N, and the reheating temperature, Trj^, are shown for some 
choices of the parameters Q, A (the values correspond to / = 1.). 
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